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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under  Program  Element  65807F. 
The  results  of  the  research  were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup 
Corporation  Company),  operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force 
Station,  Tennessee,  under  ARO  Projects  No.  V32I-A8A,  V32I-A7S,  and  V32I-A2A.  Captain 
Stanislaus  L.  Ludwig  (CF)  was  the  Air  Force  project  manager.  The  data  analysis  was 
completed  on  July  29,  1977,  and  the  manuscript  was  submitted  for  publication  on 
October  3,  1977. 
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1,0  INTRODUCTION 

The  need  for  experimental  determinations  of  the  particle  size  distribution  function 
(PSDF)  for  particulate  matter  within  the  gaseous  flow  fields  encountered  in  ground-based 
simulation  facilities  is  a continuing  requirement  to  ensure  accuracy  of  the  simulation.  The 
production  of  supersonic  and  hypersonic  flow  fields  for  aerodynamic  simulation  purposes 
normally  and  necessarily  involves  the  expansion  of  N2,  air,  or  other  gases  which  are  rarely 
free  of  particulate  matter  whether  it  be  dust,  dirt,  erosion  products,  or  condensate.  The 
presence  of  such  solid,  or  liquid  particulate  caii  result  in  both  a variety  of  deleterious 
effects  regarding  the  operation  of  the  facility  and  its  accuracy  of  simulation  and,  on  rare 
occasions,  a beneficial  effect.  Considering  the  latter  effect  first,  tlie  presence  of  either 
particulate  matter  or  condensate  within  the  flow  field  of  such  a facility  is  sometimes 
used  to  perform  so-called  unseeded  laser  velocimeter  measurements  of  the  flow  speed  and 
velocity  of  the  flow  field.  Even  for  this  example,  the  PSDF  is  of  interest  so  as  to  be 
able  to  quantitatively  predict  or  characterize  the  gas-particle  interaction  and,  specifically, 
to  determine  whether  the  particle  speed  is  an  accurate  measure  of  the  gas  flow  speed. 

The  common  effect  of  particulate  matter  in  the  flow  fields  of  ground-based  simulation 
facilities  is  one  for  which  performance  of  the  facility  is  degraded  from  its  design  conditions 
of  operation.  At  best,  the  particulate  mass  fraction  and  PSDF  may  be  such  as  to  result 
in  an  insignificant  degradation,  and  even  in  those  cases  of  observable  degradation  for  which 
the  existence  of  particles  is  known,  it  is  rarely  possible  to  determine  with  any  acceptable 
degree  of  accuracy  both  temporal  and  spatial  characteristics  of  the  PSDF  and  to  do  so 
in  a nonperturbing  manner.  In  addition  to  this  need  for  a method  to  achieve  an  accurate 
and  extensive  characterization  of  the  particulate  content  of  such  flow  fields  it  is  desirable 
that  the  same  genera!  approach  can  be  employed  to  provide  routine  monitoring  of  the 
cleanliness  of  flow  fields  whose  methods  of  production  necessarily  result  in  the  presence 
of  particulate  matter.  Examples  of  such  facilities  are  suggested  to  include  not  only  the 
more  obvious  arc-heated  flow  fields  but  also  those  which  employ  ceramic  heaters  for  which 
heater  erosion  is  known  to  occur,  facilities  which  employ  either  N2  or  air  in  either 
continuous  or  blow-down  modes  for  which  strict  gas  cleanliness  is  impossible  and,  of  course, 
all  supersonic  and  hypersonic  facilities  which  operate  within  the  test  region  either  near 
the  saturation  point  or  in  a supersaturated  condition. 

An  assessment  of  the  experimental  techniques  available  for  particle  sizing  was 
performed,  and  Mie  scattering  was  selected  as  that  technique  most  likely  to  satisfy  both 
the  then  current  and  anticipated  needs  of  the  Arnold  Engineering  Development  Center 
(AEDC)  test  facilities,  and,  consequently,  an  analytical  development  program  was  initiated 
to  determine  the  limitations  of  the  technique  and  to  ensure  the  existence  of  a theoretical 
base  which  was  adequate  for  data  analysis. 
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During  a recent  series  of  Raman  scattering  studies  of  the  flow-field  properties  of 

r t 

the  AEDC  Von  Karman  Gas  Dynamics  Facility  (VKF)  arc-heated,  hot-shot  facility  Tunnel 
F,  it  became  desirable  to  determine  the  particle  size  distribution  functions  and  the  mass 
fraction  of  particulate  matter  within  the  flow  field.  This  facility  produces  a high-enthalpy 
flow  field  of  a nominal  duration  of  130  msec,  and  both  temporal-and  spatial-resolved 
data  were  desired.  Although  all  inferential  evidence  concerning  particle  sizes  in  this  facility 
indicated  a forward  scattering  geometrical  configuration  to  be  most  desirable,  the  exigencies 
of  the  measurement  program  and  existing  tunnel  windows  dictated  the  use  of  a 
back-scattering  geometrical  configuration.  Such  measurements  were  achieved  using  two  laser 
beam  injection,  detection  channels  and  a novel  means  of  detection  of  the 
polarization-resolved  scattered  intensity  as  well  as  the  subsequent  analysis.  This  report 
presents  the  experimental  and  analytical  methods  successfully  employed,  despite  severely 
unfavorable  scattering  geometry,  in  this  initial  application  of  Mie  scattering  at  AEDC,  and 
a later  report  will  present  in  detail  the  results  of  the  measurements  and  their  implications 
as  applied  to  VKF  Tunnel  F. 

2.0  MIE  SCATTERING  AND  ITS  USE  FOR  PARTICLE  SIZE 
DISTRIBUTION  FUNCTION  DETERMINATIONS 

The  interaction  of  radiation  with  material  bodies  results  in,  among  other  processes, 
the  elastic  photon-particle  scattering  process  for  which  the  incident  wavelength  is  unchanged 
by  the  scattering  event,  and  this  process  is  depicted  in  Fig.  1 . For  particles  whose  diameters 
(D)  are  a significant  fraction  of  the  wavelength  (X)  of  the  incident  radiation  such  scattering 
is  defined  as  Mie  scattering,  which  approaches  the  Rayleigh  molecular  scattering  theory 
as  the  ratio  D/X  approaches  zero  (Ref.  1).  Restricting  consideration  to  both  scattering 
media  and  incident  fields  such  that  the  polarizability  is  a second-rank  tensor,  Mie  scattering 
theory  is  exact  for  single  scattering  events,  i.e.,  for  the  number  density  of  scatterers 
sufficiently  low  that  multiple  scattering  events  are  insignificant. 

In  addition  to  the  Mie  solution  for  spheres,  similar  exact  solutions  exist  for  scattering 
of  light  by  cylinders  (Ref.  2)  and  by  arbitrary  axisymmetric  bodies  (Ref.  3).  Further 
approximations  have  been  devised  to  treat  scattering  by  randomly  oriented  aspherical 
objects  which  do  not  fit  into  the  above  categories.  For  the  work  reported  herein  spherical 
scatterers  are  assumed. 

The  Mie  theory  assumes  the  medium  in  which  the  scatterers  are  located  is  isotropic, 
exclusive  of  the  particles  themselves.  It  is  not  necessary  to  consider  the  medium  as  being 
nondispersive  nor  nonabsorptive,  as  is  sometimes  done  in  textbooks,  and  the  Mie  scattering 
codes  we  use  are  based  on  computational  techniques  capable  of  treating  scattering  with 
complex  indices  of  refraction  for  both  particles  and  medium  (Refs.  4 and  5).  Finally, 
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it  is  assumed  that  the  number  density  of  the  molecular  medium  in  which  the  particles 
are  located  is  sufficiently  low  that  Rayleigh  scattering  by  the  surrounding  medium  is 
negligible,  and  it  is  a simple  matter  to  show  that  this  approximation  is  justified  by  all 
but  either  extraordinarily  particulate-free  gas  samples  or  gas  samples  of  extremely  high 
density.  Neither  of  these  two  cases  is  of  interest  to  this  work,  and,  consequently,  the 
medium  index  of  refraction  is  assumed  to  be  unity. 

2.1  REVIEW  OF  MIE  SCATTERING  THEORY 


The  interaction  of  collimated  radiation  of  incident  flux  I„  and  wavelength  X with 
a particle  of  diameter  D,  located  at  the  coordinate  system's  origin,  is  depicted  in  Fig. 

A A 

].  The  polarization  vector  of  the  incident  field  is  assumed  to  be  in  the  x - z plane. 


The  intensity  of  scattered  light  received  at  a point  located  a distance  d from  the 
scatterer  is  given  by  the  Mie  theory  for  each  of  the  two  polarization  states  by 


I 


1] 


.2 


—is — i,  (sin^(^>,  cos^«^) 

(kd)2 


(1) 


where  X = ttD/X.  k = 2it/X.  xt)  is  the  particle  index  of  refraction  relative  to  that  of  the 
medium,  and  6 and  0 are,  repectively,  the  polar  angle  and  azimuthal  angle  of  the 
spherical  coordinate  scheme  shown  in  Fig.  1.  Figure  2 shows  typical  variation  of  ii  and 
i2  with  X for  copper  spheres.  The  indices  1 and  2 refer,  respectively,  to  the  polarization 
states  in  which  the  scattered  light  electric  field  vector  is  perpendicular  (1 ) and  parallel  (2) 
to  the  scattering  plane.  The  angle  Q is  denoted  as  the  scattering  angle.  The  azimuthal 
factors  sin2  0 and  cos^^  are  associated,  respectively,  with  polarization  states  I 
(perpendicular)  or  2 (parallel). 


From  Eq.  (1)  the  differential  scattering  cross  section  can  be  inferred  to  be 


If  one  integrates  Eq.  2 over  47r  sr  for  each  polarization  state  and  sums  the  result,  the 
total  scattering  cross  section  (07)  is  obtained.  The  computer  program  developed  for  this 
work  gives  the  ratio  of  at  and  the  particle's  geometrical- cross  section  7r(D/2)^,  and  this 
ratio  is  defined  as  the  total  scattering  efficiency  (Qj).  The  computer  program  also  yields 
the  total  extinction  efficiency,  which  is  the  total  scattering  efficiency  plus  the  absorption 
efficiency.  The  obvious  significance  of  the  extinction  efficiency  is  that  it  determines  the 
degree  of  attenuation  of  the  radiation  as  it  passes  through  the  medium  containing  the 
scattering  particles.  If  the  optical  path  length  along  which  energy  is  lost  from  the  incident 
beam  by  scattering  and  absorption  is  U the  energy  flux  of  transmitted  light  is 


$ 


Tr 


-c 

e 


(3) 
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where  n is  the  particle  number  density  (all  sizes)  averaged  over  the  distance  L and  Oexi 
which  equals  ttQ* vt(D/2)2,  denotes  the  average  of  the  extinction  cross  section  over  the 
particle  size  distribution.  Figure  3 shows  the  variation  of  Ocxi(x)  with  particle-  size. 
Consequently,  if  the  particle  number  density  is  known,  knowledge  of  the  light  attenuation 
permits  a type  of  average  particle  size  Xext  to  be  determined.  The  precise  form  of  the 
average  is  difficult  to  determine,  but  it  is  defined  as  the  inverse  of  the  relation 


J" ^ext  ^ 


where  f(x)  is  the  particle  size  distribution  function. 


Since  the  above  described  procedure  depends- on  particle  number  density  which  is 
generally  unknown,  and  on  obtaining  accurate  absolute  intensity  measurements,  the 
procedure  is  rot  used  directly.  Instead,  a combination  of  this  procedure  with  the  angular 
scattering  procedures  is  used,  as  will  be  described  in  a following  paragraph,  and  a formal 
mathematical  deconvolution  procedure  is  applied.  Using  this  combination  enables  one  to 
infer  self-consistent  values  of  average  particle  size,  number  density  of  particles,  and 
estimates  of  particle  size  distribution  by  mole  fractions  in  each  of  as  many  as  four  size 
bins  (in  the  case  of  this  experiment)  spanning  the  size  range  of  interest.  Considering  angular 
Mie  scattering  and  assuming  that  the  solid  angle  (Af2)  is  subtended  by  a detector  located 
at  the  spherical  polar  coordinates  (d,ds,^j)  the  scattered  flux  received  by  the  detector  is 
given  by 

9,+ !i  0 ^ 

~ ~ ^ ,/*  2 ^ 

” 6^~A9/2 

where  d^2  = sin0  dd  d<p,  12  = length  of  the  scattering  volume,  and  the  ordering  of  the 
subscripts  1 and  2 and  the  parenthetical  factors  sin^0,  cos^<l>  is  to  be  maintained. 


2.2  TECHNIQUES  FOR  AVERAGE  PARTICLE  SIZE  DETERMINATIONS 

For  the  application  of  Mie  scattering  for  the  determination  of  particle  sizes  in 
aerodynamic  simulation  facilities  it  is  the  usual  case  that  the  local  total  number  density 
of  particulate  as  well  as  the  chemical  identity  of  the  scatterer  are  not  known.  For  this 
reason  and,  also,  to  eliminate  the  necessity  for  corrections  required  for  laser  power 
fluctuations,  the  approach  taken  in  this  effort  was  to  obtain  ratios  of  various  experimeritally 
determined  quantities,  and  the  Angular/ Extinction  ratio*  technique  is  an  example  of  such 
a procedure.  In  this  method  measurements  are  performed  for  the  extinction  of  the  incident 
radiation  as  well  as  for  the  power  or  photon  rate  scattered  into  detectors  located  at  various 


*This  ratio  is  related  to  a quantity  denoted  as  a phase  function  by  astrophysicists  and  radiative 
transport  theorists. 
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scattering  angles.  It  is  assumed  that  the  particles  responsible  for  the  extinction  and 
scattering  processes  are  the  same.  Using  the  transmitted  power  measurements,  one  obtains 
the  quantity  <n<iext*i>»  where  <naext*>  represents  an  average  of  the  quantity  noext® 
over  the  length  L of  the  scattering  medium.  The  additional  approximation  is  made 

where  n is  the  average  number  density  of  particles  over  the  distance  L of  the  scattering 
medium,  intervening  between  the  points  of  laser  beam  injection  and  transmission  detector. 
From  the  angular  scattering  measurement,  assuming  the  detector  to  be  located  at  scattering 
angle  0si  one  obtains  the  quantity  n£ui _2(xs,®s) 

where 


1 2 = J’a-^  {{x)  d\ 

O 

' and  n is  the  number  density  of  particles  in  the  cylindrical  scattering  volume  Vj,  with 
length  2.  It  is  assumed  that  the  ratio  nL/n2  can  be  estimated  or  measured.  Thus,  the 
experimental  quantity 


(5) 


can  be  compared  with  calculations  of  oj  ,2{x,Ss)/ocxt(x)  versus  x.  Figure  4 shows  as  a 
function  of  size  parameter  x this  calculated  cross  section  ratio  for  one  scattering  angle 
and  one  state  of  polarization.  From  Fig.  4 it  is  seen  that  such  an  experimental  measurement 
can  yield  an  estimate  of  the  range  of  particle  size  parameter  x,  or  average  particle  diameter 
Dj,  consistent  with  the  experimental  results.  It  is  obvious  that  this  method,  as  exemplified 
by  the  use  of  Fig.  4,  can  quite  conceivably  result  in  a nonunique  particle  size  parameter 
x;  e.g.,  if  the  experimental  angular/extinction  ratio  were  found  to  be,  say,  0.12,  the  range 
of  the  resulting  value  of  x is  seen  from  Fig.  4 to  be  0.5  < x < 2.4.  Consequently,  it  is 
desired  to  provide  additional  measures  for  even  though  four  plots  to  the  type  of  Fig. 
4 are  used,  corresponding  to  measurements  at  two  scattering  angles  and  in  both 
polorization  states. 


An  additional  procedure  for  providing  a separate  and  distinct  determination  of  an 
average  particle  size  parameter  involves  the  formation  of  linear  combinations  of  ratios 
of  the  angular  scattered  intensities  in  two  polarization  states.  One  such  combination, 
defined  to  be  the  degree  of  polarization  P(0s.x)  is 


P(0,.x)  = 


[ii(e,.x)]-[i^(g^.x)] 


(6) 
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where 

■■■  2 

/ i I sinOd^ 

d ~^0.'2 

S 

and  P(6x,x)  ranges  from  0 for  the  unpolarized  case  to  ±1  for  the  completely 
polarized  case. 

The  variation  of  P(0s,x)  with  size  parameter  x is  shown  in  Fig.  5 for  the  case  6^ 
- 146.5  deg.  For  a poly  disperse  sample  of  PSDF.  /(x),  one  can  define  an  average  size 
parameter  x by  the  relation 


^ (7a) 

Jt' dx 
o 

Consequently,  as  many  estimates  of  x are  provided  as  there  are  detectors  at  the  various 
scattering  angles, 

It  is  quite  obvious  that  the  resulting  value  of  x determined  by  either  the  first  method 
or  any  one  detector  using  the  second  method  is  decidedly  nonunique  as  a result  of  the 
absence  of  monotonicity  in  the  variation  of  the  dependent  function  with  x.  However, 
it  was  found  that  a combination  of  both  the  first  and  second  methods,  the  latter  employing 
two  detectors,  was  sufficient  to  reduce  the  uncertainty  in  x and  D to  acceptable  levels. 

2.3  DECONVOLUTION  OF  SIZE  DISTRIBUTION  FUNCTION 


Assuming  the  polydisperse  particulate  sample  to  be  described  by  the  PSDF,  f(x), 
the  energy  flux  ^1^2  detected  at  the  scattering  angle  0s  is  given  by 


1 ,2(^8’  “i* 


0 


(8) 


where  x^jx  is  the  largest  size  parameter  of  significance  for  the  problem  at  hand.  It  is 
obviously  of  interest  to  deconvolve  Eq.  (8)  to  extract  the  desired  PSDF,  f (x),  and  such 
a process  is  the  formal  inversion  of  Eq,  (8).  Prior  to  a discussion  of  this  inversion  procedure 
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it  is  worthwhile  to  review  the  general  nomenclature  of  integral  equations.  The  following 
equation 

b 

<^(x)  = f(x)  + A (9a) 

a 

is  known  as  either  the  Fredholm  equation  or  an  integral  equation  of  the  second  kind. 
In  Eq.  (9a)  and  the  following  Eqs.  (9)  it  is  assumed  that  t and  x are  real  continuous 
variables  and  the  functions  0(x),  f(x)  and  K(x,f)  satisfy  the  appropriate  mathematical 
requirements  of  continuity,  finiteness,  and  integrability  as  necessary  to  ensure  meaning 
for  the  equation.  In  Eq.  (9a)  K(x,^)  is  defined  to  be  the  kernel  of  the  equation,  and, 
if  K(x,£)  = K(5,x),  the  kernel  is  said  to  be  symmetric.  Further,  if  f(x)  equals  zero  one 
has  a homogeneous  integral  equation 

b 

i(x)  = A Jk(x.^)  (9b) 

a 

and,  as  in  Eq.  (9a),  X is  in  general  a complex  constant  parameter.  Finally,  the  integral 
equation  of  the  first  kind  is  defined  by 

b 

I(x)  = A rK(x,^)0(^d^  (9c) 


and,  as  in  Eqs.  (9a)  and  (9b),  a and  b are  assumed  to  be  the  real  number  limits  of  the 
integrals. 

It  is  seen  by  comparison  of  Eqs.  (8)  and  (9c)  that  the  inversion  or  deconvolution 
for  determining  f(x)  is  formally  equivalent  to  a solution  of  an  integral  equation  (Fredholm) 
of  the  first  kind.  Further,  one  can  easily  see  that  if  ratios  of  the  scattered  intensity  values 
or,  equivalently,  the  scattered  energy  flux  values  are  used  the  type  of  equation  to  be 
inverted  is  of  the  form 


111  ax 

Q 

X 

max 

/ «'i,2W„2.x)f(x)  dx 


(10) 


which  is  easily  shown  to  be  an  integral  equation  of  the  first  kind. 
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The  denominator  of  the  right-hand  side  of  Eq,  (1)  can,  of  course,  be  considered 
as  a normalization  constant.  However,  for  iterative  soiution  procedures  for  f(x)  its  value 
will  vary  as  f(x)  varies,  thereby  complicating  somewhat  the  iterative  procedure. 
Consequently,  it  is  desirable  to  select  both  a value  for  6^2  iind  a combination 
function  of  d>i  ^2(^refiX)>  say  ^ f'l’i  ,2(®rcf.x)]  such  that  ^ exhibits  an  insignificant 
variation  with  x.  If  this  is  approximately  true  one  sees  that 

X 

nt  a Y 

' = J ,2  OD 

n 

since 

X 

max 

f{x)dx  = 1. 

Therefore,  in  this  approximation,  the  value  of  ,2(®ier>x)]  assumes  the  role  of  a 
normalization  constant,  which  simplifies  the  iteration  process. 

For  the  sizing  of  particles  of  the  range  of  diameters  anticipated  for  the  present 
application,  VKF  Tunnel  F,  the  desired  configuration  of  scattering  was  estimated  to  be 
for  the  scattering  angle  range  0 < 0 < 40  deg  and  for  0tef  = 90  deg;  also  S was  anticipated 
to  be  the  sum  of  and  $2  at  flref-  However,  as  mentioned  previously,  a back-scattering 
configuration  was  required  by  other  constraints  for  this  application  and  ^ was  taken  to 
be  the  sum  of  the  scattered  flux  values  from  the  two  angular  detectors,  each  summed 
over  both  polarization  states;  i.e., 

where  0i  and  02  were  the  scattering  angles. 

The  set  of  input  data  for  the  PSDF  deconvolution  scheme  is  taken  to  be  the 
4-component  vector 


gj  = 

g2  = 4.,(02)/s 

gg  = <1)2(0])/^ 

§4  = <t>2(02^'~ 


(13) 
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The  deconvolution  scheme  is  described  in  Appendix  A and  consists  of  a matrix 
inversion  procedure  using  the  symmetric  matrix  Q^Q,  where  the  scattering  kernels 
Q are  defined  to  be 


Ij  ” 

R 

2j  = 

'Si 

|r 

O 

II 

BJR 


(14) 


where 


n = (i/N) 


(A 


j=l 


Ij 


2i 


and  N is  the  number  of  size  bins  to  be  used. 


The  quadratures  A,j  and  Bjj  are  obtained  by  numerical  integration  of  the  Mie  functions 
over  the  size  bins 


The  output  of  the  deconvolution  procedure  is  a histogram  of  mole  fractions  in  preselected 
equal  size  bins.  These  mole  fractions  are  related  to  the  actual  continuous  and  unknown 
size  distribution  function  by  Fj  = / ^F(x)dx,  and  the  size  range  xj-xj.]  defines  a size  bin. 

Xj.l 

Because  the  deconvolution  procedure  is  based  on  properties  of  the  symmetric 
scattering  kernel  matrix  M = 0^^  the  deconvolved  mole  fractions  represent  a least  squares 
solution.  Consequently,  the  requirement  that  the  number  of  size  bins  equal  the  number 
of  independent  angular  scattering  input  ratios  is  removed. 

3.0  EXPERIMENTAL  METHOD 
3.1  EXPERIMENTAL  APPARATUS 

The  experimental  arrangement  and  apparatus  employed  for  this  study  are  shown  in 
Fig.  6 and  it  is  seen  that  the  laser  beam  is  injected  diagonally  across  the  flow  field,  using 
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existing  Tunnel  F ports,  and  scattering  detectors  are  located  at  the  scattering  angles  of 
146.5  and  161.9  deg.  The  laser  source  used  was  an  Ar-Ion  Coherent  Radiation  Model 
CR-5,  source  which  was  operated  at  514.5  nm  with  an  average  output  power  of 
0.62w  as  measured  with  a Coherent  Radiation  Model  210  power  meter.  The 
laser-beam-transinitted  power  across  the  tunnel  flow  was  measured  with  a Coherent 
Radiation  Model  2 1 2 power  meter,  and  the  transmissfon  signal  was  displayed  on  a Tektronix 
Model  7633  oscilloscope  and  photographic  recording  was  used.  Since  the  level  of  visible 
fluorescence  from  the  flow  field  was  unknown  prior  to  entry,  provision  was  made  for 
electronically  modulating  the  laser  output  power  using  a signal  generator  operated  at 
nominal  1,000  Hz.  From  Fig.  6 it  is  seen  that  the  laser  beam  traversed  the  54.5-in.-diam 
tunnel  along  a window-to- window,  line-of-sight  distance  of  78.4  in.  (199  cm). 

Once  again,  because  of  our  lack  of  knowledge  of  the  flow-field  fluorescence  levels, 
spectrometers  were  employed  in  the  detection  of  the  scattered  signal  to  achieve  spectral 
isolation.  These  two  spectrometers  used  were  0.75-m  and  0.25-m  instruments,  respectively, 
and  RCA  31034A  photomultiplier  tubes  were  used  for  detection.  Data  processing  was 
accomplished  with  routine  photon  counting  instrumentation,  and  data  recording  for  the 
two  channels  was  accomplished  with  dual-beam  oscilloscopes. 

The  novel  feature  of  the  experimental  system  was  the  method  employed  to  achieve 
resolution  of  the  data  in  both  time  and  polarization.  The  collected  scattered  radiation 
was  passed  through  chopper  wheels  shown  in  Fig.  7.  and  the  modulation  frequency  was 
such  as  to  provide  a nominal  time  resolution  of  1 0 msec.  To  achieve  polarization  resolution, 
alternate  slits  in  the  chopper  were  outfitted  with  polarizer  material,  and  the  polarization 
direction  of  the  successive  slits  was  crossed.  Each  third  slit  was  blank  to  provide  an 
unpolarized  signal  as  well  as  a time  mark  reference.  Figure  7 shows  a sketch  of  this 
arrangement.  Consequently,  time-resolved  measurements  of  "ti  (6),  and  + 

<62(0)  could  be  achieved.  Finally,  a polarization  rotator  was  located  in  the  laser  injection 
optical  train  to  orient  the  input  electric  field  such  that  ^ = 45  deg,  the  importance  of  which 
is  seen  from  Eq.  (1). 

3.2  CALIBRATION  AND  OPERATION 

The  system  sensitivity  was  determined  by  performing  Rayleigh  scattering 
measurements  in  air  both  before  and  after  tunnel  firings.  Depolarization  measurements 
were  performed  to  ensure  the  absence  of  particulate  scattering  in  these  calibrations,  and 
the  known  Rayleigh  cross  sections  of  N2,  O2,  and  CO2  as  well  as  the  air  composition 
were  used  to  determine  quantitatively  the  system  sensitivity.  Further,  the  solid  angles 
subtended  by  the  detectors  at  the  scattering  volume  were  accurately  determined  to  enable 
the  conversion  from  scattered  intensity  to  power  flux. 
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TTie  firing  of  the  tunnel  provided  a trigger  source  for  the  instrumentation  channels, 
and  the  calibrated  oscilloscope  sweep  rate  as  well  as  the  chopper  modulation  frequency 
provided  time-mark  references.  Data  reduction  was  effected  in  a routine  manner,  linearity 
of  performance  of  the  overall  system  having  been  ensured. 

4.0  RESULTS 


4.1  AVERAGE  SIZE  DETERMINATION 

Using  the  previously  described  techniques,  the  measured  flux-intensity  ratios  are 

= 1.02  X 10“'  + 5% 

<I>o 

= 8.10x10“®  ± 3.6% 


and 


Oj(e2)/^o  = 2.70x10“'^  ± 6% 
= 1.68x10“'^+  4.2% 


= 0.96 


where  d\  = 146.5  deg  and  02  = 161.9  deg.  The  imprecision  of  is  estimated  to 

be  negligible,  particularly  in  view  of  the  fact  that  £n(4>xr/^o)  is  the  quantity  of  interest. 

Using  these  raw  data,  the  following  parameters  are  computed  for  use  in  the  average 
size  analysis  procedures: 


for 


Degree  of  Polarization 

P(0i)  = 0.114  ± 0.030  P(02)  = 0.232  ± 0.035 

Angular/ Extinction  Ratios 

Ri(0i)  = 0.065  ± 0.003  R2(0i)  = 0.052  ± 0.002 

Ri (02 ) = 0.116  ± 0.007  R2(02)  = 0.072  ± 0.003 


L = 48.0  in. 

cos(18.1°) 
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and 


Ri  (d, ) = 0.069  ± 0.003 
Ri(02)  = 0.123  ± 0.007 


R2(ei ) = 0.055  ± 0.002 
R2(e2)  = 0-077  ± 0.003 


for 


L = 51  .0  in. 

cosCia.l") 

The  nrst  value  of  L (representing  a lower  limit  on  the  distance  over  which  scattering 
and  absorption  occur)  is  the  nozzle  exit  diameter,  corrected  for  the  obliquity  of  the  laser 
beam.  The  second  value  of  L represents  the  obliquity-corrected  expanded  flow  dimension 
at  the  observation  station  and  was  calculated  using  a method  of  characteristics  technique. 

The  computer  simulation  of  scattering  by  various  size  distributions  of  the  form 

f{x)  = 0^ 

K! 

indicates  that  the  angular/extinction  technique  is  capable  of  average  diameter  determination 
with  perhaps  ±10%  uncertainty  if  L is  known  accurately.  In  the  present  experiment, 
knowledge  of  L is  the  limiting  factor  in  determining  the  average  diameter  accurately.  The 
range  of  angular/extinction  ratios  associated  with  the  upper  and  lower  limits  on  L has 
considerable  overlap  with  the  range  of  ratio  values  obtained  using  the  lower  limit  for 
L and  the  2-o  error  bars  of  the  experimental  data.  Consequently,  this  uncertainty  in  L 
is  within  the  imprecision  of  the  measurement  method. 

Figures  8 and  9 show  the  computed  degree  of  polarization  with  the  nominal 
experimental  values  and  error  brackets  {1  standard  deviation)  superimposed.  The  same 
type  of  graph  was  computed  for  iron  particles  and  the  results  showed  that  iron  particles 
were  not  the  major  contributor  to  the  observed  scattering.  Finally,  the  computed 
angular/extinction  ratios  are  shown  in  Figs.  1 0 through  1 3 along  with  the  measured  values 
of  these  parameters.  From  Figs.  8 through  13  it  can  be  seen  that  the  only  average  particle 
diameter  D which  is  consistent  with  the  experimental  data,  assuming  1-a  error  limits  and 
for  L = (51.0/cos  18.0,  in.  is  0.7  < D <0.9  ium.  A similar  analysis  eliminated  tungsten, 
water,  and  condensed  nitrogen  from  consideration  as  single  species  sources  of  scattering. 
Multiple  species  combinations  are  beyond  the  scope  of  this  report,  but  will  be  considered  in 
future  analyses. 
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4.2  SIZE  DISTRIBUTION  DECONVOLUTION 


Using  the  average  diameter  range  0.7  < D < 0.9  one  can  directly  determine 
the  two  bin  size  histogram,  using  the  following  relations: 


Fj  = ^~°12 


and 


P _ D-Dqi 

r 2 - 


^or°12 


where  i and  Dj  2 are,  respectively,  the  mean  diameters  corresponding  to  the  first  and 
the  second  size  bins,  and  Fi  and  F2  are,  respectively,  the  particle  mole  fractions  in  the 
first  and  second  size  bins.  Using  size  bins  encompassing  the  range  of  diameters  0 to  1.0  fxm 
and  1,0  pm  to  2.0  pm,  for  example,  we  obtain  Fj  = 0.8  and  F2  = 0.2  for  D = 0.7  pm  and  Fj 
= 0.6  and  F2  = 0.4  for  D = 0.9  pm.  Thus,  it  is  estimated  that  60  to  80  percent  of  the 
particles  have  diameters  less  than  1 .0  pm. 


More  sizing  detail  than  the  above  can  be  obtained  from  the  multibin  deconvolution 
procedure.  Computer  simulation  of  scattering  from  polydisperse  size  distributions  of  the 
form 

F(x)  = e-C* 

N! 

show  that  very  accurate  deconvolutions  are  obtained  by  the  use  of  a number  of  forward 
scattering  detector  channels  and  a 90'deg  scattering  reference  signal,  as  was  mentioned 
in  a previous  section.  Shown  in  Fig.  14  is  an  example  of  such  a deconvolution  of  simulated 
measurements  for  which  five  simulated  forward  scattering  angular  measurements  (0  to  40 
deg)  were  used  along  with  only  one  polarization  state.  Using  three  forward  and  two 
backward  angles  with  polarization  states,  we  have  recently  achieved  accurate  deconvolution 
into  twice  as  many  size  bins. 

As  stated  previously,  the  Tunnel  F scattering  geometry  is  less  than  ideal.  As  a 
consequence,  the  deconvolution  procedure  has  necessarily  been  modifled.  Deconvolution 
of  simulated  polydispersions  has  shown  that  with  these  modifications  to  the  deconvolution 
technique,  the  computed  mole  fraction  histograms  are  less  sensitive  to  the  input  data  than 
they  are  to  the  average  diameter.  Thus,  in  the  context  of  the  Tunnel  F experiment,  the 
smooth  sizing  procedure  becomes  a sophisticated  technique  for  varying  the  mole  fraction 
distribution  in  a systematic  fashion  until  the  computed  average  diameter  agrees  with  the 
average  diameter  determined  as  described  in  the  previous  section. 

Figure  15  shows  a comparison  of  the  deconvolved  Tunnel  F measurements  with  a 
simulated  polydispersion  having  a similar  average  diameter.  Although  not  shown, 
comparison  with  a similar  polydispersion  which  peaked  at  the  same  diameter  but  had 
a smaller  average  diameter  (by  a factor  of  1.5)  indicated  that  adjusting  the  smoothing 
parameter  until  correct  average  diameters  are  obtained  in  the  deconvolution  procedure 
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yields  histograms  whose  largest  mole  fractions  are  correct  to  at  least  10  to  15  percent, 
but  whose  smallest  mole  fractions  (on  the  tail  of  the  distribution)  are  unreliable.  In  fact, 
when  the  true  mole  fraction  is  less  than  one  percent,  the  computed  result  can  even  go 
negative. 

The  inputs  for  the  deconvolution  of  the  measured  light  scattering  data  are 

gi  = 0.164  g3  = 0.131 

g2  = 0.434  g4  = 0.271 

When  the  deconvolution  was  undertaken  with  zero  smoothing,  negative  mole  fractions 
were  obtained  in  two  of  the  bins  and  the  computed  average  diameter  was  an  order  of 
magnitude  too  small.  Negative  mole  fractions  did  not  disappear  until  the  smoothing 
parameter  (starting  from  zero)  reached  0.0027.  yielding  an  average  diameter  of  0.52  pm  and 
a computed  error  parameter  of  the  same  order  of  magnitude  as  the  experimental  value.  The 
mole  fraction  histogram  corresponding  to  these  values  is  given  below. 

Mole  Fraction  Bin  Diameter,  pm 


0.556 

0.008  - 

0.5 

0.332 

0.5 

0.99 

0.109 

0.9 

1.48 

0.002 

1.48  - 

1.97 

When  the  smoothing  process  is  continued  until  the  computed  average  diameter  agrees 
with  the  previously  determined  range  of  values  one  obtains  the  histogram  shown  in  Fig, 
15  and  tabulated  below  (for  smoothing  parameter  = 0.0219). 

Mole  Fraction  Bin  Diameter,  pm 


0.45 

0.008 

- 0.5 

0.31 

0.5 

- 0.99 

0.17 

0.99 

- 1.48 

0.06 

1.48 

- 1.97 

and  D = 0.67  pm. 
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The  analogous  3-bin  histogram  is: 

Mole  Fraction  Bin  Diameter,  jum 


0.58 

0.008  - 

0.66 

0.32 

0.66 

1.32 

0.10 

1.32  - 

1.97 

and  D = 0.67  pm. 

Two-bin  results  with  larger  size  range  are 


Mole  Fraction 


Bln  Diameter,  pm 


0.97 

0.03 

and  D = 0.7  pm;  and 

0.82 

0.18 


and  D = 0.9  pm. 


0.008  - 1.32 
1.32  - 2.63 


0.008  - 1.32  pm 
1.32  - 2.63  pm 


5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


Using  the  results  of  the  previous  section,  it  is  of  interest  to  obtain  an  order  of 
magnitude  estimate  of  the  average  particulate  mass  fraction  pp  which  is  defined  to  be 
the  ratio  of  the  particulate  mass  density  pp  and  the  sum  of  the  gas  density  (pg)  and 
Pp.  For  this  computation  the  particles  are  assumed  to  be  copper  spheres  of  density  p, 
= 8.89  gm/cc  and  the  flow-field  condition  is  assumed  to  be  pg  = 2.1 2 x gm/cc,  T^  = 60 
K,  and  = 2.5  x 10*  cm/s,  where  T_  and  are  the  flow-field  kinetic  temperature  and  flow 
speed,  respectively.  Using  the  transmission  relation  and  assuming  uniform  particulate 
density  across  the  flow  field  one  can  easily  show  that 

Pp/Pg  = irp^  . . £n(l/r)/6pg  . . L 

where  D is  the  mean  particle  diameter,  r is  the  transmission  factor,  op  'S  tbe  extinction 
cross  section  for  a particle  of  diameter  D,  and  L is  the  transmission  length. 

Using  the  assumed  values  one  finds 

p/p  = 66.3  (D^ /an) 

p g ° 
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and  from  Fig.  3 one  finds  that  for  D = 0.8  jum,  od  = 1.5  x 10'®  cm^. 
Consequently, 

p !p  2.3  X 10-® 

P g 

which  results  in 


“P  ■ '■p'Cp*'’,)  - (^)/('  ’ =■ 

Therefore,  the  approximate  particulate  mass  fraction  for  this  case  is  2.3  x lO"^. 


It  is  of  interest  to  recall  the  imprecision  limits  of  D : D = 0.8  ± 0.1  |um.  Translating  this 
uncertainty  in  D into  a corresponding  uncertainty  in  Op . one  finds 


showing  the  imprecision  in  ap  to  be  on  the  order  of  ±10  percent. 


Although  the  Reynolds  number  for  the  particle-gas  flow  field  cannot  be  simply 
computed,  one  can  estimate  the  likelihood  of  significant  particle  lag.  The  characteristic 
time  Tv  for  particle  lag  can  be  defined  as  (Ref.  6). 

Ty  = Dpg/18t/ 

where  v is  the  gas  viscosity  which  is  approximated  by  i/  ss  4.5  x 10"^  cgs.  From  this 
equation  for  Tv,  one  finds  for  D = 0.8  liin. 


~ 10  sec. 

Now,  a characteristic  flow  time  is  assumed  to  be  crudely  given  by  the  quantity 

where  V|h  is  the  flow  speed  at  the  nozzle  throat,  v^  is  the  free-stream  speed,  assumed 
to  be  acquired  at  the  nozzle  exit  plane,  and  AL  is  the  length  of  the  nozzle  from  throat 
to  exit.  Using  reasonable  values  for  AL,  vtj,  and  v^  one  finds 


Therefore, 


Tp  5 X 10  ® see, 

r^,  ~ 0 (100  X 


and  one  expects  that  particle  lag  effects  will  not  be  significant  since  the  flow-field  speed 
is  changing  more  slowly  than  the  speed  with  which  the  particle  can  relax  or  be 
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accommodated  to  that  change.  Consequently,  one  expects  that  the  particle  flow  speed 
Vp  “ and  the  Reynolds  number 

Re  = ■ fl  I 'to  - 

will  be  quite  small  because  of  the  near-equality  of  Vp  and  v„. 

Now,  one  can  make  use  of  this  result  to  determine  an  order  of  magnitude  estimate 
of  the  total  particle  mass  contained  within  the  flow  field  of  an  estimated  usable  run  time 
of  130  msec.  Assuming  Vp  = v_,  temporally  uniform  mass  flux,  and  that  the  circular 
cross-sectional  area  of  the  flow  field  is  1.3  x 10^  cm^,  corresponding  to  a flow  diameter 
of  — 130  cm,  one  finds  the  total  mass  (Am)  within  the  flow,  after  the  establishment 
of  the  flow,  to  be 


Am  ~ 20  gm. 

That  the  preceding  estimates  are  at  best  crude  order  of  magnitude  computations  is 
emphasized  by  the  neglect  of  the  transition  of  the  particle  flow  from  a continuum 
environment  to  a free  molecular  condition  in  the  test  section  where  the  particle  Knudsen 
number  is  estimated  to  be  on  the  order  of  100.  This  transition  is  estimated  to  occur 
between  the  nozzle  throat  and  exit  plane  and  should  be  considered  in  more  detail  for 
more  accurate  estimates. 

Other  quantitative  results  of  interest  include  the  particle  number  density.  Again,  using 
a particle  of  mean  diameter  D = 0.8  pm  for  computations  one  finds  the  mass  per  particle 
to  be  2.4  X 10"'^^  gm  and  a particle  number  density  (np)  to  be 

Hp  - 2 X 10^  cc-^ 

Further,  the  specific  volume  (Vp)  is  found  to  be 

Vp  = * 10“^  cc 

so  that  a mean  distance  between  particles  can  be  found  to  be 

(l/Hp)!/®  = Vp'-'^  ^ 3.6  X 10-2 

Finally,  using  the  well-known  criterion  for  incoherent  scattering 

2 k . sin((9/2)  ■ >>l 

one  finds  for  this  measurement 

2 k sin  (0/2)  ~10^  >>  1. 

Therefore,  the  scattering  process  is  definitely  incoherent  as  anticipated. 


21 


AEDC-TR-77-101 


In  conclusion,  this  report  has  demonstrated  both  experimental  and  analytical 
techniques  which  provide  information  concerning  the  particulate  mass  fraction  in  the  free 
stream  of  an  arc-heated  hypersonic  flow  facility.  Specifically,  an  experimental  technique 
has  been  demonstrated  which  has  the  capability  of  obtaining  time-resolved  particle  size 
distribution  functions,  and  the  time  resolution  of  this  work  was  on  the  order  of  1 0 msec, 
a value  which  can  be  substantially  decreased,  if  desired.  Further,  a novel  method  of 
measuring  in  short  times  both  states  of  polarization  of  the  scattered  radiation  has  been 
demonstrated  along  with  a method  of  analyzing  such  results  for  the  disadvantageous 
geometry  required  for  the  experimental  measurements.  Additionally,  such  results  have 
yielded  polydisperse  particle  size  distribution  functions,  mean  particle  diameter,  an  estimate 
of  the  particulate  mass  fraction  and  number  density  and  the  total  mass  of  particulate  within 
the  usable  facility  run-time.  Finally,  the  data  presented  here  represent  only  a small  portion 
of  the  flow  time  of  one  tunnel  shot.  Obviously,  the  additional  acquired  data  will  allow  a 
presentation  of  the  time-dependent  particle  size  distribution  function  as  well  as  the 
just-mentioned  mean  particle  parameters. 
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Figure  1.  Mie  scattering  experimental  configuration. 
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Figure  2.  Size  dependence  of  Mie  intensity  functions  in  geometry 
of  Fig.  1. 
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^ a.  Schematic  of  components  in  optical  path 

'O 


b.  Partial  view  of  chopper 

Figure  7.  Particle  sizing  detector  channel  configuration. 
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Degree  of  Polarization 


Figure  8.  Calculated  size  dependence  of  degree  of  polarization  with 
superimposed  range  of  measured  values,  6 = 146.5  deg. 
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Figure  9,  Calculated  size  dependence  of  degree  of  polarization  with 
superimposed  range  of  measured  values,  6 = 1 61 .9  deg. 
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Angular/Extinction  Ratio  (2) 


Figure  11.  Calculated  size  dependence  of  angular/extinction  ratio  with 
superimposed  range  of  measured  values,  0 - 146.5  deg, 
parallel  polarization. 
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Figure  12.  Calculated  size  dependence  of  angular/extinction  ratio  with 
superimposed  range  of  measured  values,  0 = 161.9  deg, 
perpendicular  polarization. 
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Figure  13.  Calculated  size  dependence  of  anguiar/extinction  ratio  with 
superimposed  range  of  measured  vaiues,  6 - 161.9  deg. 
parallel  polarization. 


AEDC-TR-77-101 


Mole  Fraction 


Diameter,  p.in 


Figure  14.  Deconvolution  of  e simulated  polydispersion. 
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Figure  15.  Comparison  of  deconvolved  Tunnel  F sizing  measurements  with 
simulated  polydispersion  deconvolution. 
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In  many  areas  of  applied  physics  it  is  impossible  to  make  direct  measurements  of 
some  desired  function  of  a quantity  x,  but,  instead,  one  must  infer  the  values  of  f(x) 
through  complementary  measurements  of  an  associated  function  of  a quantity  y.  The  two 
parameters  are  often  related  by  a Fredholm  integral  equation  of  the  first  kind. 

g(y)  = JKx)  K(y,x)dx  (A-l) 

O 

The  relation  between  particle  size  distribution  f(x)  and  scattered  light  intensity 
measurements  distributed  as  functions  of  either  wavelength  or  scattering  angle  is  described 
by  such  an  equation.  When  wavelength  is  the  argument  of  inversion  (y),  the  kernel  takes 
an  especially  simple  form  K(X,D)  =>K(D/X).  In  the  limit  of  Rayleigh-Gans  scattering,  Eq. 
(A-I)  has  been  inverted  in  closed  form  by  use  of  Mellin  transforms  (Ref.  A-l).  When 
the  more  general  case  of  Mie  scattering  is  considered,  such  an  approach  must  be  abandoned 
because  of  the  complex  dependence  of  the  Mie  functions  on  size  parameter  x = ffD/X, 
although  in  principle  the  transforms  could  still  be  obtained  numerically  be  exploiting  the 
simplified  form  of  the  kernel.  When  the  scattering  angle  is  the  argument,  there  is  no 
symmetry  in  the  dependence  of  the  kernel  on  d and  on  x,  and  the  inversion  must  be 
accomplished  by  reduction  to  quadratures  and  the  use  of  matrix  inversion  techniques. 


The  matrix  analog  of  Eq.  (A-l)  is 


- K,  F 


(A-2) 


where 


gj  = g(yi> 


f’i 


Z. 

J 


F(x)dx 


and  Ej  represents  the  combination  of  experimental  and  computational  errors  (introduced 
by  reduction  to  guadratures,  replacement  of  the  integral  of  a product  by  a product  of 
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intervals,  and  truncation  caused  by  repalcing  » ^virh  a finite  Xmax)-  Phillips  (Ref.  A-2)  has 
observed  that  direct  inversion  of  Eq.  (A-2)  neglecting  such  errors,  leads  to  unphysical 
oscillations  of  the  solution  components  Fj  — even  to  negative  values.  He  has  introduced 
a variational  technique  to  account  for  these  errors,  obtaining  a least-squares  fit  to  the 
true  solution.  Twomey  has  improved  and  generalized  Phillips'  procedure,  obtaining  solutions 
which  include  both  the  smoothest  (in  least-squares  sense)  solution  consistent  with  a given 
level  of  the  mean  square  error  and  the  solution  with  the  least  deviation  from  a trial  solution 
(Ref.  A-3). 

Dave  has  encoded  the  trial  solution  technique  and  obtained  good  results  by 
deconvolving  computer  simulated  polydispersions,  but  the  deconvolution  of  actual 
scattering  measurements  has  produced  only  fair  results  for  the  scattering  from  water 
droplets  whose  size  distribution  was  independently  determined  using  a Dessens  trap  (Ref. 
A-4).  Dave's  deconvolution  was  based  on  angular  scattering  data  obtained  and  published 
by  another  investigator,  and  there  is  some  doubt  regarding  the  quality  of  the  data. 
Glow  and  Tien  (Ref.  A-5)  have  tested  the  smoothing  technique  using  computer 
simulated  polydispersions  and  a large  number  of  simulated  angular  intensity  measurements. 
Their  work  showed  that  the  distribution  function  can  be  accurately  unfolded,  even  in 
the  presence  of  simulated  random  errors. 

Accounting  for  the  error  components  ej  and  using  the  smoothest  solution  procedure 
we  proceed  as  follows: 


Now,  let 


and  let 


setting 


gives 


IKK.)  = - 2F.  + IV, 

y = 


(A-3) 


2,  j^,u(F,)  ii'(F,)ai^  - '1  ^1-  “ 
where  y is  a lagrange  multiplier. 


(A-4) 
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Let 


k..  F. 

■'ji  1 


= Si  IKF;)  l['(Fi)5i: 


Then  it  can  be  shown  that 


^ji  = ^(j,i±2)  “ ‘'j  ^ ^ (A-5) 

The  requireinent  that  the  solution  Fj  can  be  continued  to  values  of  i less,  than  1.0  yields 

1*11  “ *’l2  = “2  = lijj,  hgg  = 5. 

Other  authors  (notably  Ref.  A-5)  have  applied  a similar  continuation  condition  at  the 
lower  right  corner  of  the  smoothing  matrix.  Such  a form  does  not  appear  in  Twomey's 
original  work,  and  it  seems  unjustifiable  to  us,  since  we  wish  to  deconvolve  the  size 
distribution  into  a comparatively  small  number  of  bins. 


Now  that  the  form  of  the  smoothing  matrix  has  been  determined,  it  is  comparatively 
simple  to  incorporate  it  into  the  solution  procedure.  Note  that 


(A-6) 


Combining  this  result  with  Eqs.  (A-4)  and  (A-3)  and  inserting  the  smoothing  matrix  yields 


+ rlk,)F,  . D\,g,  = R^,F,  (A-7) 

Therefore,  the  smoothest  least  squares  solution,  at  a given  level  of  mean  square  error 
is 


F,  - K-;  kT  g, 


The  associated  error  vector  is  obtained  from  Eq.  (A-2). 


fi  = K;.  F.  - 


with  RMS  value  t = 


and  the  relative  error  is  just 


(A-8) 


(A-9) 
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In  practice,  one  increases  7 (the  smoothing  parameter)  until  the  computed  relative  error 
exceeds  the  estimated  relative  experimental  error. 

It  is  of  some  interest  to  consider  the  significance  of  these  errors  in  order  to  understand 
the  function  of  the  smoothing  process  and  why  smoothing  is  necessary.  Twomey  has  shown 
that  unaccounted  errors  in  the  measurement  set  act  to  eliminate  the  information  content 
of  the  highest  frequency  components  of  the  Fourier  transform  of  the  kernel  (Ref.  A-6). 
In  fact,  the  dimension  of  independence  of  the  set  of  measurements  being  inverted  must 
be  reduced  by  one  for  each  eigenvalue  of  the  symmetrized  kernel  KjjjKicj  which  falls 
below  a threshold  value  related  to  the  mean  square  error  (Ref.  A-6).  In  effect,  the  kernel 
matrix  is  singular  under  these  conditions.  The  smoothing  procedure  systematically 
suppresses  the  eigenvectors  associated  with  the  occurrence  of  each  such  near-singularity. 

The  structure  of  the  symmetrized  kernel  and,  hence,  the  solution  vector  F,  changes 
drastically  when  the  smoothing  parameter  has  just  the  right  value  to  suppress  a given 
eigenvector  and  very  slowly  when  the  smoothing  parameter  has  any  other  value.  These 
critical  values  of  the  smoothing  parameter  occur  near  (but  not  necessarily  equal  to)  the 
eigenvalues  associated  with  the  eigenvectors  being  suppressed.  If  the  smoothing  matrix 
were  diagonal,  the  critical  values  of  the  smoothing  parameter  would  equal  the  associated 
eigenvalues. 


As  observed  in  the  text,  Eq.  (A-1)  is  suitable  for  inversion  of  absolute  intensity 
measurements  (as  functions  of  scattering  angle  and/or  wavelength),  but  it  must  be  modified 
when  intensity  ratios  are  used  as  the  argument  of  inversion.  Consider  a reference  kernel 
K(yo,x)  (which  can  be  a linear  combination  of  functions  Sj  a,  K(yi,x))  with  yo  chosen 
so  as  to  minimize  the  dependence  of  the  kernel  on  x.  Define  a transformed  kernel  as 


Eq.  (A- 1)  now  becomes 


where 


8(y) 


J"Q(y,K)  f(x)dx 


(A-10) 
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Normalization  is  assured  by  setting 


J f (x)  <lx 
0 

It  is  apparent  that  the  quantity  in  square  brackets  plays  the  role  of  a normalization  constant 

and  is  1 when  varies  more  slowly  with  x than  f(x). 
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NOMENCLATURE 

A Scattering  kernel  integrals  [defined  by  Eq.  (19)] 

B Scattering  kernel  integrals  [defined  by  Eq.  (19)] 

C Parameter  appearing  in  distribution  function  used  in  simulated 

deconvolution 

D Particle  diameter,  ftm 

Dqi  Mean  diameter  of  first  size  bin 

Di2  Mean  diameter  of  second  size  bin 

d Distance  from  scattering  event  to  detector,  cm 

^(1,11)  Electric  field  vectors,  perpendicular  and  parallel,  respectively,  to  scattering 

plane  (Fig.  1) 

Fs  Mole  fraction  of  particles  in  size  bin 

f(s)  Particle  size  distribution  function 

^ 4-component  vector  [defined  by  Eq.  (13)]  used  as  input  to  deconvolution 

procedure 

h_  Smoothing  matrix  defined  in  Appendix  A 

I Intensity  of  light,  w/cm^ 

i Mie  intensity  function  [defined  by  Eq.  (1)] 

K Kernel  of  integral  equation 

k Wavenumber  of  incident  light,  cm'* 

L Distance  over  which  appreciable  scattering  and  absorption  occur 

L Effective  absorption  length  used  in  angular/extinction  ratio  sizing  procedure 

£ Length  of  scattering  volume 

M Symmetric  scattering  kernel  matrix  [defined  by  Eq.  (A-3)] 
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N Parameter  appearing  in  distribution  function  used  in  simulated 

deconvolution 

n Particle  number  density 

n Particle  number  density  averaged  over  distance  L 

P(ej,x)  Degree  of  polarization  [(Eq.  (6)1 

Q Scattering  kernel  matrix  [defined  by  Eqs.  (14)  and  (15)] 

Qext  Extinction  efficiency  ratio 

Qx  Scattering  efficiency  ratio 

R Reference  denominator  for  scattering  kernel  matrix 

R Smoothed  symmetric  scattering  kernel  matrix  [defined  by  Eq.  (A-7)] 

R*  Particle  Reynolds  number  relative  to  flow 

Ri(fls)  Angular/extinction  ratio  for  "i**’"  polarization  state  evaluated  at  scattering 

angle 

V Specific  volume  occupied  by  the  average  particle 

X Particle  size  parameter,  dimensionless 

A A A 

x,y,z  Coordinate  axes  shown  in  Fig.  1 

y Dummy  variable  appearing  in  integral  inversion  discussion  in  Appendix  A 

a Particle  mass  fraction 

7 Smoothing  parameter  defined  in  Appendix  A 

AL  Distance  from  nozzle  throat  to  exit  plane 

Am  Total  particle  mass  in  flow  after  establishment  of  flow 

Ad  Range  of  scattering  angles  subtended  by  detector 

A0  Range  of  azimuthal  angles  subtended  by  detector 

e Error  vector  [defined  by  Eq.  (A-2)] 
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Subscripts 

0 

1 

2 


Particle  index  of  refraction  relative  to  medium  (in  Section  5.0  only) 
Polar  scattering  angle  (see  Fig.  1),  deg 

Wavelength  of  incident  light  in  the  scattering  medium;  also  constant  [Eqs. 
(9a)  through  (9c)]  appearing  in  discussion  of  integral  equations 

Gas  viscosity 

Dummy  variable  [Eq.  (9a)]  appearing  in  discussion  of  integral  equations 

Reference  denominator  [defined  by  Eqs.  (11)  and  (12)]  used  in 
deconvolution  procedure 

Density 

Scattering  cross  section  integrated  over  solid  angle  subtended  by  detector 
EKfferential  scattering  cross  section 
Extinction  cross  section 

Scattering  cross  section  integrated  over  47r  steradians  and  summed  over 
polarization 

Characteristic  times  appearing  in  Section  5.0 
Energy  flux  of  light  beam,  w 
Azimuthal  angles,  deg  (see  Fig.  1) 

Denotes  integration  of  enclosed  quantity  over  the  range  of  scattering  angles 
subtended  by  the  detector 

Denotes  average  of  noext^  over  the  length  L 

Incident 

Denotes  scattering  state  with  electric  field  vector  perpendicular  to  scattering 
plane 

Denotes  scattering  state  with  electric  field  vector  parallel  to  scattering  plane 
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D Denotes  a property  of  particles  of  diameter  D 

F Denotes  a property  of  the  flow 

g Elenotes  a property  of  the  gas  field 

max  Denotes  finite  value  of  x at  which  size  distribution  function  is  truncated 

p Denotes  a property  of  the  particle  field 

s Denotes  measurements  associated  with  the  set  of  angles  on  which 

a detector  is  centered 

th  Denotes  conditions  evaluated  at  the  nozzle  throat 

Tr  Denotes  transmitted  beam 

V Denotes  conditions  evaluated  at  the  mean  particle  speed 

Denotes  free-stream  conditions 

Superscripts 

Denotes  average  over  size  distribution 

^ Denotes  function  appearing  in  discussion  of  integral  equations  [Eqs.  (9a) 

through  (9c)] 

ij  Matrix  indices 
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